Probing Microquasars with TeV Neutrinos 



Amir Levinson 1 
1 School of Physics and Astronomy, Tel 
2 Department of Condensed Matter Physics, 



& Eli Waxman 2 

Aviv University, Tel Aviv 69978, Israel 
Weizmann Institute, Rehovot 76100, Israel 



o ■ 
o : 

<: 
o\ : 



(N 
> 

O 

o 



' 

a: 
i 



X 



The jets associated with Galactic micro-quasars are believed to be ejected by accreting stellar mass 
black-holes or neutron stars. We show that if the energy content of the jets in the transient sources is 
dominated by electron-proton plasma, then a several hour outburst of 1-100 TeV neutrinos produced 
by photo-meson interactions should precede the radio flares associated with major ejection events. 
Several neutrinos may be detected during a single outburst by a 1 km 2 detector, thereby providing 
a powerful probe of micro-quasars jet physics. 
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Microquasars are Galactic jet sources associated with 
some classes of X-ray binaries involving both, neutron 
stars and black hole candidates (albeit with some no- 
table differences between these two classes) Jl|-|| . During 
active states, the X-ray flux and spectrum can vary sub- 
stantially between different sub-states, with a total lumi- 
nosity that, during the so-called very high states, often 
exceeds the Eddington limit (typically a few times 10 38 
ergs s _1 in those sources). Their activity involves ejec- 
tion of jets with kinetic power that appears to constitute 
a considerable fraction of the liberated accretion energy 
(in some cases the minimum jet power inferred exceeds 
the peak X-ray luminosity @,^1), and that gives rise to 
intense radio and IR flares. Radio monitoring of some X- 
ray transients have revealed superluminal motions in cur- 
rently three of the sources , indicating that at least 
in these objects the jets are relativistic. The Lorentz fac- 
tors of the radio emitting blobs have been estimated to be 
(coincidently) T ~ 2.5 in the two superluminal sources, 
GRS 1915+105 § and GRO J1655-40 §, and some- 
what smaller (T ~ 2) in the third one, XTE J1819 @. 
Whether these are representative values or merely the re- 
sult of selection effects is unclear at present; the class of 
microquasars may contain sources with much larger T, 
rendered invisible by beaming away from us [Q . 

The temporal behavior of microquasars appears to 
be rather complex. They exhibit large amplitude vari- 
ations over a broad range of time scales and frequen- 
cies, with apparent connections between the radio, IR, 
and soft/hard X-ray fluxes [p|-^2[. The characteristics 
of the multi-waveband behavior depend on the state of 
the source, that is, whether the source is in a very high, 
soft/high or low/hard state ||. The ejection episodes 
are classified into several classes according to the bright- 
ness of synchrotron emission produced in the jet and the 
characteristic time scale of the event The duration 
of major ejection events (class A) is typically on the order 
of days, while that of less powerful flares (classes B and 
C) is correspondingly shorter (minutes to hours). The 
correlations between the X-ray and synchrotron emission 
clearly indicates a connection between the accretion pro- 



cess and the jet activity. Whether radio and IR outbursts 
represent actual ejection of blobs of plasma or, alterna- 
tively, formation of internal shocks in a quasi-steady jet is 
unclear (cf. ref. In any case, since the overall time 

scale of outbursts (minutes to days) is much longer than 
the dynamical time of the compact object (milliseconds), 
it is likely that shocks will continuously form during the 
ejection event over a range of time scales that encom- 
pass the dynamical time, owing to fluctuations in the 
parameters of the expelled wind, leading to dissipation 
of a substantial fraction of the bulk energy at relatively 
small radii. The extremely rapid variations of the X-ray 
flux often seen in these sources supports this view (e.g., 
ref [p^|). If a fraction of at least a few percent of the 
jet power is tapped for acceleration of electrons to very 
high energies, then emission of high-energy gamma-rays 
is anticipated, in addition to the observed radio and IR 
emission EGRET upper limits for some X-ray no- 

vae have been considered in this context |ll| . A recently 
discovered microquasar which appears to have a persis- 
tent radio jet jlTj], seems to coincide with an unidentified 
EGRET source having a total luminosity in excess of its 
X-ray luminosity. 

The content of jets in microquasars is yet an open is- 
sue. The synchrotron emission both in the radio and 
in the IR is consistent with near equipartition magnetic 
field, which is also implied by minimum energy consider- 
ations [Q. However, the dominant energy carrier in the 
jet is presently unknown (with the exception of the jet 
in SS433). Scenarios whereby energy extraction is as- 
sociated with spin down of a Kerr black hole favor e ± 
composition (although baryon entrainment is an issue). 
However, the pair annihilation rate inferred from the es- 
timated jet power implies electromagnetic domination on 
scales smaller than roughly 10 9 cm in the superluminal 
sources, and requires a transition from electromagnetic 
to particle dominated flow above the annihilation radius 
by some unspecified mechanism jjj . Alternatively, in sce- 
narios in which an initial rise of the x-ray flux leads to 
ejection of the inner part of the accretion disk, as widely 
claimed to be suggested by the anti-correlation between 
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the X-ray and radio flares seen during major ejection 
events (cf. ref. Q for a different interpretation), e-p 
jets are expected to be produced. A possible diagnos- 
tic of e-p jets is the presence of Doppler-shifted spec- 
tral lines, such as the Hq, line as seen in SS433. The 
detection of such lines from jets having Lorentz factor 
well in excess of unity (as is the case in the superluminal 
microquasars) may, however, be far more difficult than 
in SS433, as the lines are anticipated to be very broad 
(AA/A > 0.1). Furthermore, the conditions required to 
produce detectable flux in such sources may be far more 
extreme than in SS433 p8| . Here we propose another 
diagnostic of hadronic jets, namely emission of TeV neu- 
trinos. As shown below, for typical parameters the neu- 
trinos are produced on scales much smaller than the IR 
and radio emission and, therefore, can provide a probe of 
the innermost structure of microquasars jets. 

The picture envisaged is the following: We suppose 
that on sufficiently small scales (< 10 11 cm) a significant 
fraction of the energy liberated during ejection events 
dissipates, e.g. through the formation of internal col- 
lisionless shocks, leading to the acceleration of a non- 
thermal power law distribution of protons (and electrons) 
up to the maximum energy achievable. Detailed studies 
suggest that shock acceleration of protons may be 
highly efficient, and should produce a power law spec- 
trum, dn p (e)/de oc e~ 2 (perhaps slightly steeper in the 
case of ultra-relativistic shocks). The fraction of energy 
carried by the power law component of electrons may de- 
pend sensitively on the injection process, but generally 
thought to constitute a few percent. In the following, we 
denote by rjp the fraction of the total burst energy that is 
tapped for the acceleration of protons to nonthermal en- 
ergies. We emphasize that rj p represents essentially the 
product of the fraction of total burst energy that dis- 
sipates behind the shocks (which depends primarily on 
the duty cycle of shock formation) and the efficiency at 
which protons are accelerated to nonthermal energies by 
a single shock (that may approach 100%). Observations 
of 7-ray bursts, which are believed to arise from a simi- 
lar process of jet kinetic energy dissipation (albeit with 
different Lorentz factors and luminosities [g0[), suggest 
that a significant fraction of the jet energy is dissipated 
and converted to a power-law distribution of accelerated 
electrons |2]j| . 

As shown below, synchrotron emission by electrons 
contributes large opacity to photo-meson production over 
a range of radii that depends on the strength of the 
magnetic field. (Since the photo-meson opacity is con- 
tributed predominantly by thermal electrons, the analy- 
sis outlined below is rendered independent of the details 
of electron injection to nonthermal energies.) Between 
12% to 25% of the injected proton energy is converted, 
by virtue of the large optical depth, into muon neutrinos 
having a flat spectrum in the range between 1 and ^100 
TeV, which then freely escape the system. 



Consider a jet of kinetic power Lj — 10 38 Lj3g erg 
s , and opening angle 9 = 0.1#_i, with a corresponding 
proper energy density Uj = Lj/(ir6 2 r 2 T 2 c), propagating 
in the background of X-ray photons emitted by the ac- 
cretion disc, with luminosity L x — 10 3S L X 3s erg s -1 . In 
the scenario under consideration, significant fraction of 
the jet kinetic energy is dissipated via mildly-relativistic 
internal shocks in the jet at radii 

r > T 2 cSt = 3 x 10 7 r 2 5t_ 3 cm, (1) 

where S t = 10~ 3 <5i_3 s is the source dynamical time. 
The Thomson optical depth across the jet is given by 

t t = a T n p r9T ~ 0.05=^j» (2) 

where r — 10 8 rs cm and the proper proton number den- 
sity in the jet is n p ~ Uj/(m p c 2 ). External photons can 
therefore penetrate the jet and interact with accelerated 
protons at radii larger than roughly 10 7 cm. The pro- 
ton energy, in the jet frame, for which interaction with 
external X-ray photons is at the A-resonance is 

e pA = 0.3GeV 2 /re^ = 3 x 10 14 (re^/lkeV)- 1 eV, (3) 

where e x is the observed photon energy. 

The time available, in the jet frame, for photo-pion 
interactions is « r/Tc. Thus, the optical depth to photo- 
pion production at the A resonance, contributed by the 
external photons, is 

t P7 ~ ca peak rixw- = l ^ Lx3S (e^/lkeV)" 1 , (4) 

1 C T^g 

where cr pca k = 5 x 10~ 28 is the cross section at the 
A-resonance and n x is the number density of back- 
ground X-ray photons measured in the jet frame: n x ~ 
T(L X / e x ) / (Airr 2 c). Thus, photo-meson interactions with 
external X-ray photons at the inner-most dissipation 
radii will convert a large fraction of the accelerated pro- 
ton energy to high energy pions. 

Additional contribution to the photo-meson optical 
depth comes from synchrotron photons produced inside 
the jet by the thermal electrons. For a relativistic shock, 
the mean electron energy assuming rapid equilibration 
time is ~ 0.5(r s — l)m p c 2 , where T s is the Lorentz factor 
of the shock. For the microquasars T s — 1 ~ 1 typically. 
This then yields for the (jet frame) peak energy of the 
synchrotron spectrum, 

£7, peak — 50 - 1 J M keV, (5) 
1 r$t)-i 

where £ = 10 _1 ^_i is the equipartition parameter, de- 
fined through the relation B = yJ'&Ti^Uj. Note, that 
£7, peak reduces to infrared energies at a radius of ~ 10 13 
cm, with a corresponding light crossing time of order 
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minutes, compatible with the variability time of IR baby 
flares. 

Assuming the total energy density of synchrotron pho- 
tons to be on the order of Uj yields for the corre- 
sponding energy-dependent number density, defined as 
n syn .( e ) = e(dn syn ./de), 
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where a = 1/2 for e < e pea k, and (assuming efficient 
electron injection) a = 1 for e > e pea k- The spectrum at 
e < fpeafe results from the fast cooling of electrons. The 
ratio of electron cooling time, t syn 
the dynamical time, r/Tc, is 
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We note that at radii r < 10 9 cm, inverse-Compton scat- 
tering of electrons is suppressed since scattering is in 
the Klein-Nishsina regime for photons carrying signifi- 
cant fraction of the synchrotron energy density. At these 
radii, the fast cooling of electrons indeed implies that the 
energy density of synchrotron photons is of order Uj . At 
larger radii, the energy density of synchrotron photons 
will be somewhat suppressed due to inverse-Compton 
emission. 

A rough estimate of the energy loss rate of pro- 
tons due to photo-meson interactions gives: tz} ~ 

^pcak^s 

yn.c(Ae p /e p ), where Ae p /e p ~ 0.2 is the average 
fractional energy loss in a single collision. Equating the 
latter with the acceleration rate ~ eBc/e p , we obtain for 
the maximum proton energy, 



5 x lO^W-i) 17 



eV. 
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We emphasize that the proton energy cannot in any case 
exceed the upper limit imposed by the requirement that 
the protons be confined to the system; viz., e p < 9eBr ~ 
2 x lO^^-Ljggr -1 eV. Comparing the latter with eq. 
(0) then implies that the maximum proton energy is lim- 
ited by photo-meson losses at radii r 8 < 10 2 L 2 3S T~ 4 6z{, 
and by confinement at larger radii. 

The comoving proton energy for which interaction with 
synchrotron peak photons is at the A-resonance is 
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(9) 



Consequently, protons may be accelerated to energy ex- 
ceeding that required for photo-meson interaction with 
synchrotron peak photons for radii 



< 



x 10 
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— min^ss/r^I 7 ! 2 ) cm. (10) 



From eqs. (||) and (|T(]) it is seen that for our choice 
of parameters protons can be accelerated to energies ~ 
10 16 £_i eV in the region where photo- meson interactions 
take place. The contribution of synchrotron photons to 
the photo-meson optical depth is 
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where (3 = 1/2 for e p > e p , P cak and (3 = 1 for e p < e P]PCllk . 

Eq. ([ll]) implies that at all radii the photo-meson 
optical depth is large for protons of energy exceeding 
the threshold energy (^) for interaction with synchrotron 
peak photons. We therefore expect most of the en- 
ergy of these protons to be lost to pion production. 
Neutral pions and charged pions are produced in the 
photo-meson interaction with roughly equal probabili- 
ties, and the decay of a charged pion produces 3 neu- 
trinos, 7r + — > fi + + — > e + + v e + 17^ + v^, each carry- 
ing ~ 5% of the initial proton energy. Due to the high 
photo-meson optical depth, neutrons produced in photo- 
pion interactions will lose their energy by photo-meson 
interactions on a time scale shorter then the neutron life 
time. This implies similar production rates of both neg- 
atively and positively charged pions and, hence, similar 
production rates of and 9^. 

High energy pions and muons may loose significant 
fraction of their energy by inverse-Compton scattering 
prior to decaying. Inverse-Compton scattering of syn- 
chrotron photons by high energy pions is in the Klein- 
Nishina regime for 
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0-4 , r rr f~! 1/2 TeV. (12) 



Comparing this equation with Eq. (H), and noting that 
the energy of a pion produced by photo-pion interaction 
is ~ 20% of the initial proton energy, we find that pions 
produced by protons of energy exceeding the threshold 
energy satisfy > e^^. Thus, the pion inverse-Compton 
scattering energy loss time is 



tic 



e-7r^7,pcak 



Sca T n syn . (m e c 2 ) 2 



= 3 x 10-^-i 



ITeV 



(13) 



and the ratio of the pion decay time to energy loss time 
is 



tic 

7"dccay 



(14) 



We therefore expect pions to decay prior to significant 
energy loss. However, muons, the life time of which is 
w 100 times longer, may lose significant fraction of their 
energy before decaying. In the following we therefore 
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conservatively assume that a single high energy (or 
Vfj) is produced in a single photo-pion interaction of a 
proton (or neutron), corresponding to conversion of 1/8 
of the energy lost to pion production to muon neutrinos. 
In terms of r] p , the fraction of jet energy injected as a 
power law distribution of protons, the flux at Earth of 
and can be expressed as, 



AnD 2 

= 0.5xl0-\_ 1 T- 1 5 3 D2*L j38 ergs^cm" 2 , (15) 



where 8 is the Doppler factor, D = 10 22 Z?22 cm is the dis- 
tance to the source, and rj p = 0.1?7p._i. The factor 1/2 
is due to the fact that for a flat power-law proton dis- 
tribution, and for the characteristic parameters invoked, 
roughly half the energy of the power-law component is 
carried by protons which lose most of their energy to 
pion-production [see Eqs. (pjll)]. 

Due to the high photo-meson optical depth at energies 
ak, we expect a neutrino spectrum similar to the 
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proton spectrum, dn v /de v cx e~ 2 , above ~ 1 TeV. The 
probability that a muon neutrino will produce a high- 
energy muon in a terrestrial detector is P V{1 w 1.3 x 
l(r 6 £^ TeV , with (3 = 2 for E UtTeV < 1 and /3 = 1 for 
-E"y,TeV > 1- Thus, for a flat neutrino spectrum above 
1 TeV, the muon flux at the detector is « {Pq/Eq)^ , 
where P Q /E = 1.3 x lO^TeV -1 . The number of events 
detected in a burst of energy E = 10 43 i?43 erg is therefore 



A M ~ 0.2r )p ,- 1 r- 1 S 3 D2iE 43 (A/lkm 2 ), 



(16) 



where A is the effective detector area. For a persistent 
source (like SS433), the corresponding rate of neutrino 
induced muon detection is 



0.2r/ p ,„ir~Vi^2 Lj,38(A/lkni )day~ 



(17) 



For sources directed along our sight line r -1 <5 3 ~ 8r 2 
may exceed 100. Thus, if the fraction r\ v exceeds a 
few percent, several neutrinos can be detected during 
a typical outburst even from a source at 10 kpc. The 
typical angular resolution of the planned neutrino tele- 
scopes at TeV energies (e.g. Q) should be 9 ~ 1 degree. 
The atmospheric neutrino background flux is § v ,bkg ~ 
10~ 7 e m ?Tcv/ cm2ssr ' implying a number of detected 
background events N^kg ~ 3 x lO _2 (0/deg) 2 idaykm -2 
per angular resolution element over a burst duration 
ltday day. The neutrino signals above ~ 10 TeV from 
a typical micro-quasar outburst should therefore be eas- 
ily detected above the background. 

The duration of the neutrino burst should be of the 
order of the blob's ejection time. It should precede the 
associated radio outburst, or the emergence of a new su- 
per luminal component, that originate from larger scales 
(~ 10 15 cm), by several hours. 



As an example, let us consider the 1994 March 19 out- 
burst observed in GRS 1915+105 §. This source is at a 
distance of about 12 kpc, or Z? 2 2 = 3.5. From the proper 
motions measured with the VLA, a speed of 0.92c and 
angle to the line of sight of 9 ~ 70° are inferred for the 
ejecta, corresponding to a Doppler factor 8 ~ 0.58. A 
conservative estimate of the total energy released during 
this event yields E i3 = 20 From eq. @ we then 

obtain, ~ 0.03^-1 for a 1 km 2 detector. At >10 
TeV this is above the background, but would require an 
average of ~ 30 outbursts for detection. We note that if 
a similar event were to occur in a source that is located 
closer to Earth, at a distance of say 3 kpc as in the case 
of the superluminal source GRO J1655-40 (which has a 
similar Doppler factor), or from a jet oriented at a smaller 
angle, then a few or even a single outburst may produce 
several muon events in such a detector. This suggests 
that even sources that are beamed away from us might 
be detectable during particularly strong outbursts. 

As a second example consider the source SS433. This 
source, at a distance of 3 kpc, exhibits a steady jet that 
moves with a speed of ~ 0.3c. The presence of H a lines 
indicates a baryonic content. A conservative estimate 
of the jet's kinetic power yields Lj > 10 39 ergs s _1 [ pi[ . 
Thus of order 10 3 77 Pj _i events per year per km 2 are antic- 
ipated. This rate is consistent with the MACRO upper 
limit [ p5[ . Emission of TeV neutrinos from this source 
has been proposed earlier |26[ , although the mechanism 
invoked in this ref. is pp collisions. 

Positive detection of TeV neutrinos from microquasar 
jets would imply a baryonic content (althogh non- 
detection does necessarily imply a pair dominated jet), 
and would have important consequences for the mecha- 
nisms responsible for the ejection and confinement of jets. 
It would also provide important information concerning 
dissipation and particle acceleration in shocks. Contem- 
poraneous detections of gamma-rays may even enable us 
to infer the relative efficiencies at which electrons and 
protons are accelerated in shocks by comparing the to- 
tal energies emitted as neutrinos and gamma-rays during 
outbursts, thereby resolving a long standing issue. 
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